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Abstract 

Lunar laser ranging (LLR) is used to conduct high-precision measurements of ranges between 
an observatory on Earth and a laser retro-reflector on the lunar surface. Over the years, LLR has 
benefited from a number of improvements both in observing technology and data modeling, which led 
to the current accuracy of post-fit residuals of ^ 2 cm. Today LLR is a primary technique to study 
the dynamics of the Earth-Moon system and is especially important for gravitational physics, geodesy 

■ and studies of the lunar interior. When the gravitational physics is concerned, LLR is used to perform 
, high- accuracy tests of the equivalence principle, to search for a time- variation in the gravitational 

constant, and to test predictions of various alternative theories of gravity. The gravitational physics 
parameters cause both secular and periodic effects on the lunar orbit that are detectable with the 
present day LLR; in addition, the accuracy of their determination benefits from the 35 years of the 
LLR data span. On the geodesy front, LLR contributes to the determination of Earth orientation 
parameters, such as nutation, precession (including relativistic precession), polar motion, and UT1, 
i.e. especially to the long-term variation of these effects. LLR contributes to the realization of both 
the terrestrial and selenocentric reference frames. The realization of a dynamically defined inert ial 
reference frame, in contrast to the kinematically realized frame of VLBI, offers new possibilities for 
mutual cross-checking and confirmation. Finally, LLR also investigates the processes related to the 
Moon's interior dynamics. Here, we review the LLR technique focusing on its impact on relativity 
and give an outlook to further applications, e.g. in geodesy. We present results of our dedicated 

■ studies to investigate the sensitivity of LLR data with respect to the relativistic quantities; we also 
present the computed corresponding spectra indicating the typical periods related to the relativistic 

VO . effects. We discuss the current observational situation and the level of LLR modeling implemented 

to date. We emphasis the need for the modeling improvement for the near future LLR opportunities. 
\ We also address improvements needed to fully utilize the scientific potential of LLR. 
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1 Introduction 



Being one of the first space geodetic techniques, lunar laser ranging (LLR) has routinely provided 
observations for more than 35 years. The LLR data are collected as normal points, i.e. the combination 
of lunar returns obtained over a short time span of 10 to 20 minutes. Out of ^ 10 19 photons sen t 
per pulse by the transmitter, less than 1 is statistically detected at the receiver Wil liams et al.l (fl996); 
this is because of the combination of several factors, namely energy loss (i.e. 1/R A law), atmospherical 
extinction and geometric reasons (rather small telescope apertures and reflector areas). Moreover, the 
detection of real lunar returns is rather difficult as dedicated data filtering (spatially, temporally and 
spectrally) is required. These conditions are the main reason, why only a few observatories worldwide 
are capable of laser ranging to the Moon. 

Observations began shortly after the first Apollo 11 manned mission to the Moon in 1969 which 
deployed a passive retro-reflector on its surface. Two American and two French-built reflector arrays 
(transported by Soviet spacecraft) followed until 1973. 1 Most observations are taken to the largest 

■"^One of the reflector-arrays (of the Soviet Luna 17 mission, see also Fig. |8j has not been tracked operationally. The 
reason could be that the coordinates are not known well enough or that the reflectors are covered by dust or the transport 
cap. 
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Figure 1: (a) Left: Lunar observations per year, 1970 - 2005. (b) Data distribution as a function of the 
synodic angle D. 



reflector array, that of the Apollo 15 mission. Over the years more than 16,000 LLR measurements have 
by now been made of the distance between Earth observatories and lunar reflectors. Most LLR data 
have been collected by the Observatoire de la Cote d'Azur (OCA, France), the McDonald observatory 
(Texas, USA) and - until 1990 - Haleakala (Hawaii, USA). The new data are still coming, but today only 
the first two stations operate regularly. Understanding unexpected and small effects is very difficult with 
only one or two operating stations, because possible instrumental systematics of the ranging system can 
not be separated from real scientific effects reliably. In order to further increase the impact of LLR in 
Relativity and Earth sciences more stations, with a wide geographical distribution, are needed. Therefore 
the Italian colleagues have set up a new site in Matera which has provided first LLR data quite recently. 
A new site with lunar capability is currently being built at the Apache Point Obse rvatory, New Mexico, 
USA. This station, cal led APOLLO, is designed for a mm-level accuracy ranging ((Murphy et all EoOCt 
Willi anis et allEo04bl) . However, to fully exploit the available LLR potential, a few more sites capable 
of tracking the Moon are needed, especially at diverse locations including the Southern hemisphere. 

Fig.^i shows the number of LLR normal points per year since 1970. As shown there and in Fig.^, 
the range data have not been accumulated uniformly; substantial variations in data density exist as 
a function of synodic angle D, these phase angles are represented by 36 bins of 10 degree width. In 
Fig. 03, data gaps are seen near new Moon (0 and 360 degrees) and full Moon (180 degrees) phases, and 
asymmetry about quarter Moon (90 and 270 degrees) phases also is exhibited. The former properties 
of this data distribution are a consequence of operational restrictions, such as difficulties to operate 
near the bright sun in daylight (i.e. new Moon) or of high background solar illumination noise (i.e. 
full Moon). Note also asymmetry about quarter Moon phases. The uneven distribution with respect 
to the lunar sidereal angle shown in Fig. [2Jl represents the increased difficulty of making observations 
from northern hemisphere observatories to the Moon when it is located over the southern hemisphere. 
Here, the situation will only be improve if a southern observatory (e.g. in Australia) will start to track 
the Moon. It might be possible that new missions to the Moon could be helpful in this respect; the 
deployment of active laser transponders could allow satellite laser ranging systems to participate in LLR. 

While measurement precision for all model parameters benefit from the ever-increasing improvement 
in precision of individual range measurements (which now is at the few cm level, see also Fig. 
some parameters of scientific interest, such as time variation of Newton's coupling parameter G/G or 
precession rate of lunar perigee, particularly benefit from the long time period (35 years and growing) 
of range measurements. 

In the 1970s LLR was an early space technique for determining Earth orientation parameters (EOP). 
Today LLR still competes with other space geodetic techniques, and because of large improvements in 
ranging precision (30 cm in 1969 to ~ 2 cm today), it now serves as one of the strongest tools in the 
solar system for testing general relativity. Moreover, parameters such as the station coordinates and 
velocities contributed to the International Terrestrial Reference Frame ITRF2000. EOP quantities were 
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Figure 2: (a) Data distribution as a function of the sidereal angle 5, where the Moon is south of equator 
from 0° to 180°. (b) Weighted residuals (observed-computed Earth-Moon distance) annually averaged. 

used in combined solutions of the International Earth Rotation and Reference Systems Service IERS 
(a = 0.5 mas). 

2 LLR Model and Relativity 

The existing LLR model at IfE has been developed to compute the LLR observables — the roundtrip 
travel times of laser pulses between stations on the Earth and passive reflectors on the Moon (see e.g. 
Miiller et al. 1996, Miiller and Nordtvedt 1998 or Miiller 2000, 2001, Miiller and Tesmer 2002, Williams 
et al. 2005b and the references therein). The model is fully relativistic and is complete up to first 
post-Newtonian (1/c 2 ) level; it uses the Einstein's general theory of relativity - the standard theory of 
gravity. The modeling of the relativistic parts is much more challenging than, e.g., in SLR, because the 
relativistic corrections increase the farther the distance becomes. The basic observation equation reads 

d = c- = |r em - ration + r reflect or| + c At (1) 

where d is the station-reflector distance, c the speed of light, r the pulse travel time, r em the vector 
connecting the geocenter and the selenocenter, r sta tion the geocentric position vector of the observatory, 
and r re fl ect0 r the selenocentric position vector of the reflector arrays. Ar describes corrections of the 
travel time caused by atmospheric effects, but also (relativistic) transformations into the right time 
system as well as the light time equation (Shapiro effect). In order to apply Eq. QJ, all vectors have to 
be transformed in one common reference frame (in our case the inertial frame) which requires consistent 
relativistic transformations, so-called pseudo-Lorentz transformations. The Earth-Moon vector r em can 
only be obtained by numerical integration of the corresponding equation of motion (here in simplified 
version): 

I*em — 3 F em LxlVlg I o 3 I ^newtonian ~r Drelativistic- \^) 

r em \l r s — r em| r s J 

r em is the relative acceleration vector between Earth and Moon, GM e -^ m is the Earth-Moon mass times 
the gravitational constant, r s the geocentric position vector of the Sun, r s the Earth-Sun distance, M s 
the solar mass, b newt0 nian comprises all further Newtonian terms like the effect of the other planets or 
the gravitational fields of Earth and Moon, and b re i at i v istic indicates all 'relativistic' terms, i.e. those 
entering the Einstein-Infeld-Hofmann (EIH) equations. Corresponding relativistic equations are applied 
to describe the rotational motion of the Moon. The rotation angles are then applied to rotate the 
selenocentric reflector coordinates of Eq. (1) into the inertial frame. For the transformation of the 
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geocentric station coordinates, the well known rotation matrices (precession, nutation, GAST, etc.) are 
used, see IERS 2003. 

The modeling of the 'Newtonian' parts has been set up according to IERS Conventions (IERS 2003), 
but it is restricted to the 1 cm level. The weights are based upon the accuracy estimates for the normal 
points as provided by the observatories. Based upon this model, two groups of parameters (170 in total) 
are determined by a weighted least-squares fit of the observations. The first group includes the so-called 
'Newtonian' parameters such as 

• geocentric coordinates of three Earth-based LLR stations and their velocities; 

• a set of EOPs (luni-solar precession constant, nutation coefficients of the 18.6 years period, Earth's 
rotation UT0 and variation of latitude by polar motion); 

• selenocentric coordinates of four retro-reflectors; 

• rotation of the Moon at one initial epoch (physical librations); 

• orbit (position and velocity) of the Moon at this epoch; 

• orbit of the Earth-Moon system about the Sun at one epoch; 

• mass of the Earth-Moon system times the gravitational constant; 

• the lowest mass multipole moments of the Moon; 

• lunar Love number and a rotational energy dissipation parameter; 

• lag angle indicating the lunar tidal acceleration responsible for the increase of the Earth-Moon 
distance (about 3.8 cm/yr), the increase in the lunar orbit period and the slowdown of Earth's 
angular velocity. 

The second group of parameters is used to perform LLR tests of viable modifications of the general 
theory of relativity. The general theory of relativity is not expected to be perfect, because Einstein's 
theory and quantum mechanics are fundamentally incompatible. Therefore, it is important in physics 
to find out at what level of accuracy it fails. Relativistic parameters to be determined by LLR analyses 
are (values for general relativity are given in parentheses): 

1. Strong equivalence principle (EP) parameter 77, which for metric theories is 77 = 4/3 — 3 — 7 (= 0). 
In LLR a violation of the equivalence principle would show up as a displacement of the lunar orbit 
along the direction to the Sun. LLR is the dominant test of the strong equivalence principle, i.e. 
for self-gravitating bodies. 

2. Space-curvature parameter 7 (= 1) and non-linearity parameter (3 (= 1). LLR also has the capa- 
bility of determining the PPN parameters (3 and 7 directly from the point-mass orbit perturbations 
(i.e. as described by the EIH equations), but, e.g., (3 may be derived much better by combining 
the EP parameter 77 with an independent determination of 7 (see below). 

3. Time variation of the gravitational coupling parameter G/G (= yr _1 ). This is the second most 
important gravitational physics result that LLR provides. Einstein's theory does not predict a 
changing G, but some other theories do. So it is important to measure this as well as possible. 
The sensitivity improves like the square of the data span. 

4. Geodetic de Sitter precession f^s of the lunar orbit (p± 1.92 "/cy). LLR has also provided the only 
accurate determination of the geodetic precession. The dedicated space mission GP-B will provide 
an improved accuracy, if that mission is successfully completed. 

5. Coupling constant a (= 0) of Yukawa potential for the Earth-Moon distance which corresponds to 
a test of Newton's inverse square law. 

6. ol\ (= 0) and (\2 (= 0) which parametrize 'preferred frame' effects in metric gravity. 
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7. Combination of parameters ("1 — Co — 1 (= 0) derived in the Mansouri and Sexl (1977) formalism 
indicating a violation of special relativity (there: Lorentz contraction parameter £i = 1/2, time 
dilation parameter Co = —1/2). 

8. EP-violating coupling of normal matter to 'dark matter' at the galactic center. It is a similar test 
to item 1 above, but now different periods are affected (mainly the sidereal month). 

9. A further application is the detection of the Sun's J2 10 -7 ) from LLR data (independent to 
other methods as solar seismology), which affects the anomalous perihelion shift of Mercury, one 
of the classical tests of relativity. The present uncertainty (« 10 -6 ) is larger than the expected 
value. The parameter J 2 is not further discussed in this paper and will be addressed in more detail 
in an upcoming study. 

The determination of the relativistic parameters indicated above can be accomplished either by 
modifying the equations of motion (i.e. parametrizing present terms or adding new ones) or by deriving 
analytical expressions for their effect on the Earth-Moon distance. In the first case the needed partial 
derivatives can be computed by numerical differentiation, in the second case by direct derivation of the 
analytical terms. 

Many relativistic effects produce a sequence of periodic perturbations of the Earth-Moon range 

n 

Ar EM = [ A t cos (uiAt + fy) + BiAt sm(uiAt + $;) . (3) 
i=i 

Ai and Bi, uji, and fa are the amplitudes, frequencies, and phases, respectively, of the various perturba- 
tions. Some sample periods of perturbations important for the measurement of various parameters are 
given in Table \? 



Table 1: Typical periods of some relativistic quantities, taken from Miiller et al. (1999). 



Parameter 


Typical Periods 


77 


synodic (29 d 12 /l 44 m 2.9 s ) 


G/G 


secular + emerging periodic 




sidereal, annual, sidereal-2- annual, 
anomalistic (27 d 13H8 m 33.2 s ) ± annual, synodic 




2-sidereal, 2 -sidereal- anomalistic, nodal (6798 d ) 


Ci - Co - 1 


annual (365. 25 d ) 


^galactic 


sidereal (27 d 7 /l 43 m 11.5 s ) 



3 Sensitivity Study 

As indicated in Eqs. (l)-(2), LLR is affected in various ways and at various levels by relativity. Relativity 
enters the equation of motion, i.e. the orbit and the rotation of the Moon. More precisely, the Earth- 
Moon system behaves according to relativity. But also the light propagation and the transformations 
between the reference and time frames which are used have to be modeled in agreement with general 
relativity. The lunar measurements contain the summed signal of all effects in one, so that the separation 
of the individual effects is a big challenge. To better understand what the individual contributions of 
the different relativistic effects are, sensitivity studies have been carried out, as 

ArP - S JjHL Ap . (4) 



era 



5p 



2 Note: the designations should not be used as formulae for the computation of the corresponding periods, e.g. the 
period 'sidereal-2-annual' has to be calculated as 1/(1/27. 32 d — 2/365. 25 d ) « 32.13 d . In addition, 'secular + emerging 
periodic' means the changing orbital frequencies induced by G/G are starting to become better signals than the secular 
rate of change of the Earth-Moon range in LLR. 
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Figure 3: (a) Sensitivity of LLR with respect to G/G assuming AG/G = 8 • 10 -13 yr _1 . (b) Sensitivity 
of LLR with respect to space curvature 7, non-linearity couplings (3 and geodetic precession using their 
present LLR accuracy (see Table 2) as perturbation value. 




Days since December 27, 1969 

Figure 4: (a) Sensitivity of LLR with respect to Yukawa interaction parameter a, equivalence principle 
violation 77 and time- variable gravitational constant G/G using their present LLR accuracy (see Table 
2) as perturbation value, (b) Sensitivity of LLR with respect to preferred-frame effects ql\ and 0L2 using 
their present LLR accuracy (see Table 2) as perturbation value. 

Ar v ern is the perturbation of the Earth-Moon distance caused by a parameter p which is one of the 
relativistic parameters described in Section 2. Sr em /Sp is the partial derivative of the Earth-Moon 
distance with respect to p; it is obtained by numerical differentiation. Ap is a small value indicating the 
variation in p, here we have used the present realistic error as derived from LLR analyses (see Table 2). 
As an example, Fig. Oh represents the sensitivity of the Earth-Moon distance with respect to a possible 
temporal variation of the gravitational constant in the order of 8 • 10 -13 yr _1 , the present accuracy of 
that parameter. It seems as if perturbations of up to 9 meters are still caused, but this range (compared 
to the ranging accuracy at the cm level) can not fully be exploited, because the lunar tidal acceleration 
perturbation is similar. Fig. 03, Fig. Eh and Fig. QJd show the results of corresponding investigations for 
all relativistic parameters which were investigated during the present study (i.e. without the parameters 
discussed in items 7-9 of the previous Section). The perturbations vary between 5 cm (ai) and 25 m (/?) 
which indicates that the former parameter is determined quite well from LLR data, as the sensitivity 
values are close to the observation accuracy, whereas the latter is only poorly determined because of its 
high correlation with the Newtonian orbit perturbations. Nevertheless, the continuation of LLR over a 
longer time span will help to further de-correlate the various parameters. 
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Figure 5: (a) Power spectrum of a possible equivalence principle violation assuming A(mc / raj) « 10 13 
(b) Power spectrum of the effect of G/G in the Earth-Moon distance assuming AG/G = 8 • 10 -13 yr _1 . 




Frequency [1 /d] Frequency [1 t d] 



Figure 6: (a) Power spectrum of an additional (deviating from Einstein's value) geodetic precession 
assuming Agpm = 10 -2 . (b) Power spectrum of a possible Yukawa term using Aa = 2 • 10 -11 . 

To better understand those couplings, corresponding power spectra have been computed. The largest 
periods for the EP-parameter are shown in Fig. EH and for G/G in Fig. 03. Obviously many periods 
are affected simultaneously, because the perturbations, even if caused by a single beat period only (e.g. 
the synodic month for rj parameter 3 ), change the whole lunar orbit (and rotation) and therefore excite 
further frequencies. For comparisons also the spectra of the geodetic precession gpm and the Yukawa 
coupling parameter a are indicated (Fig. EH and Fig. Et>). Again a different combination of periods is 
visible. As before mainly the monthly (e.g. sidereal, synodic, anomalistic), half- monthly, etc., periods 
dominate, but longer periods (low frequencies) are also present, e.g. annual, 3 years or combinations of 
the monthly and annual frequencies. Similar pictures are obtained when considering the preferred- frame 
parameters a\ and 0L2 (Fig. EH and Fig. Ub)- A huge variety of periods show up again, but they differ 
partly from each other (note, e.g., the very low- frequency contributions). The spectra of 7 and /3 (not 
shown here) are very similar to the geodetic precession spectrum. Although there are big similarities 
between the various spectra, the typical properties of each can be used to identify and separate the 
different effects and to determine corresponding parameters. 

3 Note, here and throughout the paper, the relation mc/mj = 1 + r](U/Mc 2 ) has been used equivalently, where the 
second term describes the self energy of a body, cf. Williams et al. 2005b. 
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Figure 7: (a) Power spectrum of a possible preferred- frame effect by ql\ assuming Aai =9-10 5 . (b) 
Power spectrum of a possible preferred- frame effect by a 2 assuming A«2 = 2.5 • 10 -5 . 

4 Results 

The global adjustment of the model by least-squares-fit procedures gives improved values for the es- 
timated parameters and their formal standard errors, while consideration of parameter correlations 
obtained from the covariance analysis and of model limitations lead to more 'realistic' errors. Incom- 
pletely modeled solid Earth tides, ocean loading or geocenter motion, and uncertainties in values of fixed 
model parameters have to be considered in those estimations. See Table 2 for the determined values for 
the relativistic quantities and their realistic errors. 

The EP-parameter 77 (= (6 ± 7) • 10 -4 ) benefits most from highest accuracy over a sufficient long 
time span (e.g. one year) and a good data coverage over the synodic month, as far as possible. Also any 
observations reducing the asymmetry about the quarter Moon phases (compare Fig. H]d) would improve 
the fit of 77. The improvement of the EP parameter was not so big since 1999, as the LLR RMS residuals 
increased a little bit in the past years, see Fig. EJd. The reason for that increase is not completely 
understood and has to be investigated further. 

In combination with the recent value of the space-curvature parameter 7cassini 
(7 — 1 = (2.1 ± 2.3) • 10 -5 ) derived from Doppler measurements to the Cassini spacecraft (Bertotti et al. 
2003), the non-linearity parameter (3 can be determined by applying the relationship tj — 4/3 — 3 — 7Cassini- 
One obtains f3 — 1 = (1.5 ± 1.8) • 10 -4 (note that using the EP test to determine parameters 77 and /3 
assumes that there is no composition- induced EP violation). The LLR result for the space-curvature 
parameter 7 as determined from the EIH equations is less accurate than the results derived from other 
measurements, because its effect is very similar to the Newtonian orbit perturbations. 

For the temporal variation of the gravitational constant, G/G = (6 ± 8) • 10 -13 has been obtained, 
where the formal standard deviation has been scaled by a factor 3 to yield the given value. This parameter 
benefits most from the long time span of LLR data and has experienced the biggest improvement over 
the past years (cf. Miiller et al. 1999). 

For the estimation of the de Sitter precession of the lunar orbit, a Coriolis-like term is added to the 
equation of motions, which adds the precession effect as predicted by Einstein for a second time. This 
term is scaled by a parameter called gpm which has to give if Einstein is correct, gpm = is obtained 
with an accuracy of about 1 percent. 

The Yukawa coupling parameter a has been determined by adding corresponding perturbation terms 
to the equations of motion, where the partials were computed by numerical differentiation. The recently 
determined value shows small deviations from the expected value; it will be further investigated in future. 

The preferred-frame parameters ql\ and 0L2 can either be determined by extending the equations of 
motion or by adding analytical terms to the Earth-Moon distance. In both cases quite similar results 
are obtained (see Miiller et al. 1996, 1999). Recent determinations are given in Table 2. 
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The Mansouri-Sexl parameters Co and fi as well as the quantity indicating a possible EP-violating 
coupling with dark matter were not investigated during our present studies; the values given in Table 2 
are taken from Miiller et al. (1999). 

A further question to be considered in more detail in future is the right combination of the parameters. 
That means, shall all relativist ic parameters be estimated together in one global adjustment, or each one 
alone (together with the parameters of the standard solution)? We carried out several tests considering 
the correlation of the relativistic quantities with each other, but also with the 'classical' ones, e.g. with 
the orbital parameters or site velocities (Koch 2005). It is too early to make a final decision. On the 
one hand 'over-estimation' of an effect has to be avoided, on the other hand over-constraining by fixing 
too many parameters should also be avoided. 

Final results for all relativistic parameters obtained from the IfE analysis are shown in Table 2 
(see also Miiller et al. 2005). The realistic errors are comparable with those obtained in other recent 
investigations, e.g. at JPL (see Williams et al. 1996, 2004a, 2004b, 2005b). 

Table 2: Determined values for the relativistic quantities and their realistic errors. 



Parameter 


Results 


Equivalence Principle parameter n 


(6 ± 7) • 10- 4 


Metric parameter 7 — 1 


(4 ±5) • 10- 3 


Metric parameter [5—1: direct measurement 
and from rj = 4(3 - 3 - 7c aS sini 


(-2 ±4) • 10- 3 
(1.5±1.8)-1(T 4 


Time varying gravitational constant G/G [yr _1 ] 


(6±8) • 10- 13 


Differential geodetic precession Qqp - ^deSit ["/ c y] 


(6 ±10) -10- 3 


Yukawa coupling constant a (for A = 4 • 10 5 km) 


(3 ±2) • 10- 11 


'Preferred frame' parameter a\ 


(-7 ±9) • 10- 5 


'Preferred frame' parameter 0L2 


(1.8 ±2.5) • 10- 5 


Special relativistic parameters ("1 — Co — 1 


(-5± 12) • 10- 5 


Influence of dark matter 5^ ga iactic [cm/s 2 ] 


(4±4)-10- M 



5 Further Applications 

In addition to the relativistic phenomena discussed above, more effects related to lunar physics, geo- 
sciences, and geodesy can be investigated. The following items are of special interest (see also Miiller et 
al. 2005), as they touch recent activities in the afore- mentioned disciplines: 

1. Celestial reference frame: A dynamical realization of the International Celestial Reference System 
(ICRS) by the lunar orbit is obtained (a = 0.001") from LLR data. This can be compared and 
analyzed with respect to the kinematical ICRS from Very Long Baseline Interferometry (VLBI). 
Here, the very good long-term stability of the lunar orbit is of great advantage. 

2. Terrestrial reference frame: The results for the station coordinates and velocities, which are esti- 
mated simultaneously in the standard solution, contribute to the realization of the international 
terrestrial reference frame, e.g. to the last one, the ITRF2000. 

3. Earth rotation: LLR contributes, among others, to the determination of long-term nutation pa- 
rameters, where again the stable, highly accurate orbit and the lack of non-conservative forces 
from air-drag or solar radiation pressure (which affect satellite orbits substantially) is very conve- 
nient. Additionally UT0 and VOL values are computed, which stabilize the combined EOP series, 
especially in the 1970s when less data from other space geodetic techniques were available. The 
precession rate is another example in this respect. 

4. Relativity: As discussed in the previous sections, the dedicated investigation of Einstein's theory 
of relativity is of major interest. With an improved accuracy of the LLR measurements and the 
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modeling (see next Section) the investigation of further effects (e.g. the Lense-Thirring precession) 
or those of alternative theories might become possible. 



5. Lunar physics: By the determination of the libration angles of the Moon, LLR gives access to 
underlying processes affecting lunar rotation (e.g. Moon's core, dissipation), cf. Williams et al. 
(2005a). A better distribution of the retro-reflectors on the Moon (see Fig. |H| would be very 
helpful. 

6. Selenocentric reference frame: The determination of a selenocentric reference frame, the combina- 
tion with high-resolution images and the establishment of a better geodetic network on the Moon 
is a further big item, which then allows accurate lunar mapping. 

7. Earth-Moon dynamics: The mass of the Earth-Moon system, the lunar tidal acceleration, possible 
geocenter variations and related processes as well as further effects can be investigated in detail. 

8. Time scales: The lunar orbit can also be considered as a long-term stable clock so that LLR can 
be used for the independent realization of time scales, which can then be compared or combined 
with other determinations. 

Those features shall be addressed in the future, when more and better LLR data are available and 
the analysis models have been improved to the mm level, see next Section. 



6 Model and Observation Refinement 

To exploit the full available potential of LLR, the theoretical models as well as the measurements require 
optimization. Using the 3.5 m telescope at the APOLLO site in New Mexico, USA, a mm-level ranging 
becomes possible. To allow an order of magnitude gain in determination of various quantities in the 
complete LLR solution, the current models have to be up-dated according to the IERS conventions 2003, 
and made compatible with the IAU 2000 resolutions. This requires, e.g., to better model 

• higher degrees of the gravity fields of Earth and Moon and their couplings; 

• the effect of the asteroids (up to 1000); 

• relativistically consistent torques in the rotational equations of the Moon; 

• relativistic spin-orbit couplings; 




Figure 8: Distribution of retro-reflectors on the lunar surface. 
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• torques caused by other planets like Jupiter; 

• the lunar tidal acceleration with more periods (diurnal and semi-diurnal); 

• ocean and atmospheric loading by updating the corresponding subroutines; 

• nutation using the recommended IAU model; 

• the tidal deformation of Earth and Moon; 

• Moon's interior (e.g. solid inner core) and its coupling to the Earth-Moon dynamics. 

Besides modeling, the overall LLR processing shall be optimized. The best strategy for the data fitting 
procedure needs to be explored for (highly) correlated parameters. 

Finally LLR should be prepared for a renaissance of lunar missions where transponders or new retro- 
reflectors may be deployed on the surface of the Moon which would enable many pure SLR stations to 
observe the Moon. NASA is planning to return to the Moon by 2008 with Lunar Reconnaissance Orbiter 
(LRO), and later with robotic landers, and then with astronauts until 2018. The primary focus of these 
planned missions will be lunar exploration and preparation for trips to Mars, but they will also provide 
opportunities for science, particularly if new reflectors are placed at more widely separated locations than 
the present configuration (see Fig.JHJ). New installations on the Moon would give stronger determinations 
of lunar rotation and tides. New reflectors on the Moon would provide additional accurate surface 
positions for cartographic control (Williams et al. 2005b), would also aid navigation of surface vehicles 
or spacecraft at the Moon, and they also would contribute significantly to research in fundamental and 
gravitational physics, LLR-derived ephemeris and lunar rotation. Moreover in the case of co-location 
of microwave transponders, the connection to the VLBI system may become possible which will open a 
wide range of further activities such as frame ties. 

7 Conclusions 

LLR has become a technique for measuring a variety of relativistic gravity parameters with unsurpassed 
precision. Sensitivity studies have been performed to estimate the order of magnitude of relativistic 
effects on lunar ranges and the potential capability to better determine certain relativistic features. 
Spectral analyses showed the typical frequencies related to each effect, indicating as well, how (highly) 
correlated parameters might be separated. No violations of general relativity are found in our investiga- 
tions. Both the weak and strong forms of the EP are verified, while strong empirical limitations on any 
inverse square law violation, time variation of G, and preferred frame effects are also obtained. 

LLR continues as an active program, and it can remain as one of the most important tools for testing 
Einstein's general relativity theory of gravitation if appropriate observations strategies are adopted and 
if the basic LLR model is further extended and improved down to the millimeter level of accuracy. The 
deployment of transponders on the Moon would considerably improve the performance for lunar ranging 
applications. Lunar science, fundamental physics, control networks for surface mapping, and navigation 
would benefit. Demonstration of active devices would prepare the way for very accurate ranging to Mars 
and other solar system bodies. 
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